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Abstract. The analysis of anisotropics in the cosmic microwave background (CMB) has become 
an extremely valuable tool for cosmology. We even have hopes that planned CMB anisotropy 
experiments may revolutionize cosmology. Together with determinations of the CMB spectrum, 
they represent the first cosmological precision measurements. This is illustrated in the talk by 
Anthony Lasenby. The value of CMB anisotropics lies to a big part in the simplicity of the theoret- 
ical analysis. Fluctuations in the CMB can be determined almost fully within linear cosmological 
perturbations theory and are not severely influenced by complicated nonlinear physics. 

In this contribution the different physical processes causing or influencing anisotropics in the 
CMB are discussed. The geometry perturbations at and after last scattering, the acoustic oscilla- 
tions in the baryon-photon-plasma prior to recombination, and the diffusion damping during the 
process of recombination. 

The perturbations due to the fluctuating gravitational field, the so called Sachs-Wolfe con- 
tribution, is described in a very general form using the Weyl tensor of the perturbed geometry. 



1 Introduction 



The formation of cosmological structure in the universe, inhomogeneities in the matter dis- 
tribution like quasars at redshifts up to z ~ 5, galaxies, clusters, super clusters, voids and 
walls, is an outstanding basically unsolved problem within the standard model of cosmol- 
ogy. We assume, that the observed inhomogeneities formed from small initial fluctuations 
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by gravitational clustering. 



At first sight it seems obvious that small density enhancements can grow sufficiently 
rapidly by gravitational instability. But global expansion of the universe and radiation 
pressure counteract gravity, so that, e.g., in the case of a radiation dominated, expanding 
universe no density inhomogeneities can grow significantly. Even in a universe dominated 
by pressure-less matter, cosmic dust, growth of density perturbations is strongly reduced by 
the expansion of the universe. 

Furthermore, we know that the universe was extremely homogeneous and isotropic at 
early times. This follows from the isotropy of the 3K Cosmic Microwave Background (CMB), 
which represents a relic of the plasma of baryons, electrons and radiation at times before 
protons and electrons combined to neutral hydrogen. After a long series of upper bounds, 
measurements with the DMR instrument aboard the COsmic Background Explorer satellite 
(COBE) have finally established anisotropics in this radiation Q at the level of 

/ _ 10 -io Qn angular scales 7 ° < q < go" . 



Such an angle independent spectrum of fluctuations on large angular scales is called 
Harrison Zel'dovich spectrum |J. It is defined by yielding constant mass fluctuations on 
horizon scales at all time, i.e., if Zjy(t) denotes the expansion scale at time t, 

((AM/M) 2 (X = l H )) = const. , independent of time. 

The COBE result, the observed spectrum and amplitude of fluctuations, strongly support 
the gravitational instability picture. 

Presently, there exist two main classes of models which predict a Harrison-Zel'dovich 
spectrum of primordial fluctuations: In the first class, quantum fluctuations expand to super 
Hubble scales during a period of inflationary expansion in the very early universe and 'freeze 
in' as classical fluctuations in energy density and geometry (see also the contribution by 
V. Mukhanov) . In the second class, a phase transition in the early universe, at a temperature 
of about 10 16 GeV leads to topological defects which induce perturbations in the geometry 
and in the matter content of the universe fij. Both classes of models are in basic agreement 
with the COBE findings, but differ in their prediction of anisotropies on smaller angular 
scales. 

On smaller angular scales the observational situation is at present somewhat confusing 
and contradictory |5|, but many anisotropies have been measured with a maximum of 
about AT/T « (3 ± 2) x 10~ 5 at angular scale 9 « (1 ± 0.5)°. There is justified hope, 
that the experiments planned and under way will improve this situation within the next few 
years (see contribution by A. Lasenby) In Fig. 1, the experimental situation as of spring '96 
is presented. 

In this paper we outline a formal derivation of general formulas which can be used to 
calculate the CMB anisotropies in a given cosmological model. Since we have the chance 
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Figure 1: The corresponding quadrupole amplitude Q/i a t is shown versus the corresponding spher- 
ical harmonic index i. The amplitude Q fi a t {£) corresponds roughly to the temperature fluctuation 
on the angular scale 9 ~ tt/£. The solid line indicates the predictions from a standard cold dark 
matter model. (Figure taken from ref. [5]). 

to address a community of relativists, we make full use of the relativistic formulation of the 
problem. In Section 2 we derive Liouville's equation for massless particles in a perturbed 
Friedmann universe. In Section 3 we discuss the effects of non- relativistic Compton scattering 
prior to decoupling. This fixes the initial conditions for the solution to the Liouville equation 
and leads to a simple approximation of the effect of collisional damping. In the next Section 
we illustrate our results with a few simple examples. Finally, we summarize our conclusions. 



Notation: We denote conformal time by t. Greek indices run from to 3, Latin indices 
run from 1 to 3. The metric signature is chosen ( — h ++). The Friedmann metric is thus 
given by ds 2 = a 2 {t){—dt 2 +^ i jdx' L dx^) 1 where 7 denotes the metric of a 3-space with constant 
curvature K. Three dimensional vectors are denoted by bold face symbols. 
We set H = c = k Bo itzmann = 1 throughout. 



2 The Liouville equation for massless particles 
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2.1 Generalities 



Collision-less particles are described by their one particle distribution function which lives 
on the seven dimensional phase space 

V m = {(x,p) e TM\g{x)(p,p) = -m 2 } . 

Here A4 denotes the spacetime manifold and TM. its tangent space. The fact that collision- 
less particles move on geodesies translates to the Liouville equation for the one particle 
distribution function, /. The Liouville equation reads |jj 

X g (f) = . (2.1) 

In a tetrad basis (e At )^ =0 of Ai, the vector field X g on V m is given by (see, e.g., ]7)) 

X 9 = (p%-o; i »^) , (2.2) 



where uo u ^ are the connection 1-forms of (Ai, g) in the basis e M , and we have chosen the basis 

d_ 

'dp 



(e M )J =0 and (— ) \ =l on TV m , P = P% 



We now show that for massless particles and conformally related metrics, 

_ 2~ 

(X g f)(x,p) =0 is equivalent to (Xgf)(x,ap) = . (2.3) 
This is easily seen if we write X g in a coordinate basis: 

with ^ 

The variables 6 M are the components of the momentum p with respect to the coordinate basis: 

P = P% = & M <9 M • 

If (e M ) is a tetrad with respect to g, then = ae^ is a tetrad basis for g. Therefore, the 
coordinates of of ap = ap^e^ = a 2 p fJl e^ = a 2 b fJl d fJi , with respect to the basis <9 M on (A4,g) are 
given by a 2 b^. In the coordinate basis thus our statement Eq. ( |2.3j ) follows, if we can show 
that 

{X g f){x^a 2 b l ) = iff (X g f)(x^b l ) = (2.4) 

Setting v = ap = v^e^ = w^d^, we have = ap^ and = a 2 b^. Using p 2 = 0, we obtain 
the following relation for the Christoffel symbols of g and g: 

r a/3 b a bP = r al3 b a bP + ^b a u . 



For this step it is crucial that the particles are massless! For massive particles the statement 
is of course not true. Inserting this result into the Liouville equation we find 



a 2 X g f = w^d,f\ b - 2^g) - f><V^ , (2.5) 
where d^f\b denotes the derivative of / w.r.t. x M at constant (b l ). Using 



d,f\ b = dj\ w + 2^b^, 
we see, that the braces in Eq. ( |2.5|) just correspond to d^f\ w . Therefore, 

a 2 XJ(x,p) = w ^f\ w -r a( ,w a w^^- = Xgf(x,ap) , 

which proves our claim. This statement is just a precise way of expressing conformal invari- 
ance of massless particles. 



2.2 Free, massless particles in a perturbed Friedmann universe 

We now apply this general framework to the case of a perturbed Friedmann universe. For 
simplicity, we restrict our analysis to the case K =, i.e., Q — 1. The metric of a perturbed 
Friedmann universe with density parameter Q — 1 is given by ds 2 = g liv dx fl dx v with 

flV = a 2 (r/^ + h^) = a 2 g^ v , (2.6) 

where (77^) = diag(—, +, +, +) is the flat Minkowski metric and (h^) is a small perturbation, 

From Eq. fl2.3|) , we conclude that the Liouville equation in a perturbed Friedmann uni- 
verse is equivalent to the Liouville equation in perturbed Minkowski space, 

(X- g f)(x,v)=0, (2.7) 

with v = v^e^ = ap M e M .[] 



We now want to derive a linear perturbation equation for Eq. fl2.7|) . If e M is a tetrad 
in Minkowski space, = + \h v ^t v is a tetrad w.r.t the perturbed geometry g. For 
(x, v^e^) G Po, thus, (x, f M e M ) G Po- Here P denotes the zero mass one particle phase space 
in Minkowski space and Po is the phase space with respect to g, perturbed Minkowski space. 
We define the perturbation, P, of the distribution function by 

f{x, v"e M ) = f{x, u%) + Fix, u%) . (2.8) 



■"■Note that also Friedmann universes with non vanishing spatial curvature, K ^ 0, are conformally flat 
and thus this procedure can also be applied for K 7^ 0. Of course, in this case the conformal factor a 2 is 
no longer just the scale factor but depends on position. A coordinate transformation which transforms the 
metric of K 7^ Friedmann universes into a conformally flat form can be found, e.g., in J^]. 
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Liouville's equation for / then leads to a perturbation equation for F. We choose the natural 
tetrad 

e =3 - -h y 8 
with the corresponding basis of 1-forms 

Inserting this into the first structure equation, dO^ = — v A fix 1 ', one finds 

Ufiu = —-^{h^\, u —h u \^)6 x . 

Using the background Liouville equation, namely that / is only a function of v — ap, we 
obtain the perturbation equation 

(d t + n l di)F = -~[{h i0 - h 00 ,i )n % + (toy - h oj ,i K nJ ]^; > 

where we have set t>j = t>nj, with t> 2 = X^iX^) 2 ; i- e -> n gives the momentum direction of 
the particle. Let us parameterize the perturbations of the metric by 

{kflu) = ( Bi m L b tJ + 2H tJ ) ' (2 - 9) 
with Hj = 0. Inserting this above we obtain 

(d t + n%)F = -[H L + (A,i + l -R i )rf + (Hij - ^A,>W>£ . (2.10) 

From Eq. (|2.10|) we see that the perturbation in the distribution function in each spectral 
band is proportional to v-f-. This shows once more that gravity is achromatic. We thus do 
not loose any information if we integrate this equation over photon energies. We define 

7i r q _ 

m = \ Fv dv . 

p r a A J 

Am is the fractional perturbation of the brightness l, 

l = a~ A / fv 3 dv . 



Setting t(n, x) = t(T(n, x)), one obtains that i = (7r/60)T 4 (n, x). Hence, m corresponds to 
the fractional perturbation in the temperature, 

T(n,x) = f (1 +m(n,x)) . (2.11) 

Another derivation of Eq. fl2.11| ) is given in |TD| . According to Eq. fl2.1U| ) , the v dependence 
of F is of the form vj£. Using now 

4tt J ^v A dv = -A J fv 3 dvdn = -Ap r a A , (2.12) 



we find 

df 

F(x^,n\v) = -m(x M ,n*)f— . 

dv 

This shows that m is indeed the quantity which is measured in a CMB anisotropy experiment, 
where the spectral information is used to verify that the spectrum of perturbations is the 
derivative of a blackbody spectrum. Of course, in a real experiment located at a fixed 
position in the Universe, the monopole and dipole contributions to m cannot be measured. 
They cannot be distinguished from a background component and from a dipole due to our 
peculiar motion w.r.t. the CMB radiation. 

Multiplying Eq. (|2.10| ) with v 3 and integrating over v, we obtain the equation of motion 
for m 

d t m + n%m = H L + + ^B l )n l + {H i3 - -B uj )n i n j . (2.13) 

It is well known that the equation of motion for photons only couples to the Weyl part 
of the curvature (null geodesies are conformally invariant). However, the r.h.s. of Eq. ( j2. 13 ) 



is given by first derivatives of the metric only which could at best represent integrals of the 
Weyl tensor. To obtain a local, non integral equation, we thus rewrite Eq. fl2.13|) in terms 



of V 2 m. It turns out, that the most suitable variable is however not V 2 m but x, which is 
defined by 

x = V 2 m - (V 2 H L - ) - ^(V 2 ^; - 3dV 4 ,K , 

where cry = -\^B Uj +B jH ) + -S^Bf + H i5 . 

Note that \ and V 2 m only differ by the monopole contribution, V 2 Hl — (l/2)H lJ ,y , and 
the dipole term, (l/2)(V 2 5j — ?>d^aij)n l . The higher multipoles of \ and V 2 m agree. An 
observer at fixed position and time cannot distinguish a monopole contribution from an 
isotropic background and a dipole contribution from a peculiar motion. Only the higher 
multipoles, / > 2 contain information about temperature anisotropies. For a fixed observer 
therefore, we can identify V~ 2 x with 5T/T. 

In terms of metric perturbations, the electric and magnetic part of the Weyl tensor are 
given by (see, e.g. 0, [Toj ) 

Sij = i[A lJ (A-i7 L )-a lJ -V 2 ^-^^ + ^+4 n ] ( 2 - 14 ) 

with A y = didj - (l/3)5y V 2 . 

Explicitly working out (d t +n l di)x using Eq. fl2.13p , yields after some algebra the equation 
of motion for x' 

(d t + n%) X = Zrt&Bij + T^nUmdiBij = S(t, x, n) , (2.16) 
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where e^u is the totally antisymmetric tensor in three dimensions with €123 = 1. The spatial 
indices in this equation are raised and lowered with 5y and thus index positions are irrelevant. 
Double indices are summed over irrespective of their positions. 

Eq. (|2.16|) is the main result of this paper. We now discuss it, rewrite it in integral form 
and specify initial conditions for adiabatic scalar perturbations with or without seeds. 



In Eq. ( |2.16|) the contribution from the electric part of the Weyl tensor is a divergence, 



and therefore does not contain tensor perturbations. On the other hand, scalar perturbations 
do not induce a magnetic gravitational field. The second contribution to the source term 
in Eq. Q2.16 ) thus represents a combination of vector and tensor perturbations. If vector 



perturbations are negligible (like, e.g., in models where initial fluctuations are generated 
during an epoch of inflation), the two terms on the r.h.s of Eq. ( [OBJ ) yield thus a split into 
scalar and tensor perturbations which is local. 



Since the Weyl tensor of Friedmann Lemaitre universes vanishes, the r.h.s. of Eq. ( [2.16| ) 



is manifestly gauge invariant (this is the so called Stewart-Walker lemma |12|). Hence also 



the variable x 1S g au g e invariant. Another proof of the gauge invariance of x-, discussing the 



behavior of F under infinitesimal coordinate transformations is presented in flip 
The general solution of Eq. ( j2.16| ) is given by 



X(t,x,n)= t S(f,x + (t' -t)n,n)dt' + x(U,x + (U - t)n,n) , (2.17) 

Jti 

where S is the source term on the r.h.s. of Eq. (|2.16|) . 

In Appendix A we derive the relations between the geometric source term S and the 
energy momentum tensor in a perturbed Friedmann universe. 



3 The collision term 

In order for Eq. (|2.17|) to provide a useful solution, we need to determine the correct initial 
conditions, x(tdec), at the moment of decoupling of matter and radiation. Before recom- 
bination, photons, electrons and baryons form a tightly coupled plasma, and thus x can 
not develop higher moments in n. The main collision process is non-relativistic Compton 
scattering of electrons and photons. The only non vanishing moments in the distribution 
function before decoupling are the zeroth, i.e., the energy density, and the first, the energy 
flow. We therefore set 

X(tdec) = V 2 (~!>< r) (W) - n • VW(W)) , (3.1) 

where 

D { ;\t dec ) = V- 2 J X (tdec)d^j (3.2) 



5 -^—-AH L + 2V- 2 {H >) and 
P 
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V j( - r \t dec ) = -V- 2 (— /x(Wy^) (3.3) 

= -T (ry /(^W) + - 2v- a (^ y ) . 

DW and V {r) are gauge invariant density and velocity perturbation variables || |10[| . 

In the tight coupling or fluid limit, the initial conditions can also be obtained from the 
collision term. Setting M. = V~ 2 x one finds the following expression for the collision integral 

[0, 

C[M] = aa T n e [-D^ -M + n- V {b) + -n^M' 11 ] . (3.4) 
4 y 2 

The last term is due to the anisotropy of the cross section for non-relativistic Compton 
scattering, with 

M n = — [ (rf n i _ -gV)Mdn . 
8n J 3 

M. is a gauge invariant perturbation variable for the distribution function of photons. 
denotes the baryon velocity field, ot and n e are the Thomson cross section and the free 
electron density respectively. To make contact with other literature, we note that Ai = G+$, 
where G is the perturbation variable describing the CMB anisotropies defined in [D| and 



$ denotes a Bardeen potential (see Section 4). Since Ai and G differ only by a monopole 
term, they give rise to the same spectrum of temperature anisotropies for £ > 1. Ai satisfies 
the Boltzmann equation 

(d t + n l di)M = V- 2 S + C[M] , (3.5) 

where S is the gravitational source term given in Eq. ( |2.16| ). In the tight coupling limit, 
t T = (aaTn e )~ l t, we may, to lowest order in (tx/t), just set the square bracket on the 
right hand side of Eq. Q3.4p equal to zero. Together with Eq. ( |3.3| ) this yields 

y(6) _ y(r) 

Neglecting gravitational effects, the right hand side of Boltzmann's equation then leads to 

£>(r) = ^ v . = 4^(6) j (3 6) 

where the last equal sign is due to baryon number conservation. In other words, photons and 
baryons are adiabatically coupled. Expanding Eq. ( |3.5| ) one order higher in tx, one obtains 
Silk damping |fT4|, the damping of radiation perturbations due to imperfect coupling. 



Let us estimate this damping by neglecting gravitational effects and the time dependence 
of the coefficients in the Boltzmann equation ( |3.5|) since we are interested in time scales 
t T t. We can then look for solutions of the form 

V {b) oc M oc exp(i(kx - ujt)) . 
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We also neglect the angular dependence of the collision term. Solving Eq. ( |3.5| ) for Ai, we 
then find 

(l/4)D^ + ik ■ nV® 

M = V 1 . 9 , ; — . 3.7 

1 - it T (u -k-n) v ; 

The collisions also induce a drag force in the equation of motion of the baryons which is 
given by 

Fi = f C[M]n t dn = - ikV®) . 

With this force, the baryon equation of motion becomes 

kujV {b) + i(a/a)kV {b) = ikm - F / p h . 

To lowest order in tx/t and Ht, this leads to the following correction to the adiabatic 
condition V {b) = V {r) : 

t T ukV® = ^tikV® - V«) , (3.8) 
3p& 

From Eq. (|3.6|) we obtain the relation k • = —(?>/ A)ujD^ to lowest order. Using this 
approximation, we find, after multiplying Eq. ( |3.8j ) with k, 



with R = 3pb/p r - The densities pb and p r denote the baryon and radiation densities respec- 
tively. Inserting this result in Eq. ( |3.7|) leads to 

■y _|_ 3flLd/k 

M = . + Y tTUjR u r ^ r) /4 , (3.10) 

1 — it T (uj — hp:) y 

where we have set pi — k • n/k. From this result, which is valid on time scales shorter 
than the expansion time (length scales smaller than the horizon), we can derive a dispersion 
relation oj{k). In lowest order utx we obtain 

uj = ujq — «7 with (3-11) 

u Q = and 7 = kh T 5 . 3.12 

^3(1 + R) Q{R + 1) 2 

At recombination R ~ 0.1 so that 7 ~ 2k 2 tx/15. 

We have thus found that, due to diffusion damping, the photon perturbations thus un- 
dergo an exponential decay which can be approximated by 

\M.\ oc exp (-2k 2 t T t/15) , on scales t > l/k > t T . (3.13) 



In general, the temporal evolution of radiation perturbations can be split into three 
regimes: Before recombination, t td ec the evolution of photons can be determined in the 
fluid limit. After recombination, the free Liouville equation is valid. Only during recombi- 
nation the full Boltzmann equation has to be considered, but also there collisional damping 
can be reasonable well approximated by an exponential damping envelope ||15|| , which is a 
somewhat sophisticated version of (|3.13|) . 
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4 Example: Adiabatic scalar perturbations 



We now want to discuss Eq. ( |2.16| ) with initial conditions given by Eq. ( |3.1| ) in some examples. 

Perturbations are called 'scalar' if all 3 dimensional tensors (tensors w.r.t their spatial 
components on hyper-surfaces of constant time) can be obtained as derivatives of scalar 
potentials. 

Scalar perturbations of the geometry can be described by two gauge invariant variables, 
the Bardeen potentials |16| $ and The variable \1/ is the relativistic analog of the Newto- 
nian potential. In the Newtonian limit, — $ = \]/= the Newtonian gravitational potential. In 
the relativistic situation, $ is better interpreted as the perturbation in the scalar curvature 
on the hyper-surfaces of constant time ||17|| . In terms of the Bardeen potentials, the electric 
and magnetic components of the Weyl tensor are given by [Tl[ 

1 

2' 



Sij = -A y ($ - *) , 2% = 0, (4.1) 



where Ay denotes the traceless part of the second derivative, Ay = didj — |<5jjV 2 . The 
Liouville equation, fl2.16p then reduces to 

(d t + n%)M = n l di(<$> - V) . (4.2) 

With the initial conditions given in Eq. ( |3.1| ) we find the solution 

S -I(t ,x ,n) =M{t ,Xo,n) = []d^ +n%V® -$}{t dec , x dec )- f° ($-*)(*, x{t))dt , 
J- 4 Jt dec 

(4.3) 

where x dec = x — (t Q — td ec ) n an d correspondingly x(t) = x — (t — t)n. 

We now want to replace the fluid variables, and V^ b \ wherever possible, by perturba- 
tions in the geometry. To this goal, let us first consider the general situation, when one part 
of the geometry perturbation is due to perturbations in the cosmic matter components and 
another part is due to some type of seeds, which do not contribute to the background energy 
and pressure. The Bardeen potentials can then be split into contributions from matter and 
seeds: 

$ = $ m + $ s , * = ^ m + ^ s . (4.4) 

To proceed further, we must assume a relation between the perturbations in the total energy 
density and energy flow, D g and V, and the corresponding perturbations in the photon 
component. The most natural assumption here is that perturbations are adiabatic, i.e., that 

D^/(l + Wr ) = D g /(l + w) and V {b) = V {r) = V , 

where w = p/p denotes the enthalpy, i.e. w r = 1/3. For w r ^ w this condition can only 
be maintained on super-horizon scales or for tightly coupled fluids. For decoupled fluid 
components, the different equations of state lead to a violation of this initial condition on 
sub-horizon scales. 
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In order to use the perturbed Einstein equations to replace D g and V by geometric 
perturbations we define yet another density perturbation variable, 

D = D g + 3(1 + w)-V -3(1 + w)$ and 

(X 

D (r) = D (r) + A ^y(r) _ 4$ _ 

y a 

The matter perturbations D and V determine the matter part of the Bardeen potentials via 
the perturbed Einstein equations (see, e.g. [fLCR ). The following relation between $ m and D 
can also be obtained using Eqs. ( |4.1| ) and (A16) in the absence of seeds. 

D = -\(~) V 2 $ m ~ (H) 2 $ m and 

-* m -* m = - - i+™v. 

a 2 Va/ 

The term D rsp. £)^, is much smaller than the Bardeen potentials on super-horizon scales 
and it starts to dominate on sub-horizon scales, ^ > 1. For this term therefore, the 
adiabatic relation is not useful and we should not replace by 3 ^ W ^ D. The same 

holds for diV^ which is of the order of kt<& m . However, (d/a)V A( - r - ) is of the same order 
of magnitude as the Bardeen potentials and thus mainly relevant on super horizon scales. 
There the adiabatic condition makes sense and we may replace (a/a)V by its expression 
in terms geometric perturbations. Keeping only and diV^ in terms of photon fluid 
variables, Eq. (|4.3|) becomes 

5 -^(x ,t ,n) = [g, + l±ftE ro+ — | (-) 1 $ m + ^W+nWW](x fc W) 
1 3 + 6w 6(1 + w ) \a/ 4 

(4.5) 

This is the most general result for adiabatic scalar perturbations in the photon temper- 
ature. It contains geometric perturbations, acoustic oscillations prior to recombination and 
the Doppler term. Silk damping, which is relevant on very small angular scales (see the 
contribution by 0) is neglected, i.e., we assume 'instantaneous recombination'. Eq. ( |4.5|) is 
valid for all types of matter models, with or without cosmological constant and/or spatial 
curvature (we just assumed that the latter is negligible at the last scattering surface, which 
is clearly required by observational constraints). The first two terms in the square bracket 
are usually called the ordinary Sachs-Wolfe contribution. The integral is the 'integrated 
Sachs- Wolfe effect'. The third and fourth term in the square bracket describe the acoustic 
Doppler oscillations respectively. On super horizon scales, kt <C 1, they can be neglected. 

To make contact with the formula usually found in textbooks, we finally constrain our- 
selves to a universe dominated by cold dark matter (CDM), i.e., w = without any seed 
perturbations. In this case \I/ S = $ s = and it is easy to show that ^ = — $ and that 
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$ = ^ = (see, e.g., ]T0|). Our results then simplifies on super-horizon scales, kt <C 1, to 
the well-known relation of Sachs and Wolfe [|ll| 

= -V(x -t n,t dec ) . (4.6) 
sw 6 




5 Conclusions 



We have derived all the basic ingredients to determine the temperature fluctuations in the 
CMB. Since the fluctuations are so small, they can be calculated fully within linear cos- 
mological perturbation theory. Note however that density perturbations along the line of 
sight to the last scattering surface might be large, and thus the Bardeen potentials inside 
the Sachs Wolfe integral might have to be calculated within non-linear Newtonian gravity. 
But the Bardeen potentials themselves remain small (as long as the photons never come 
close to black holes) such that Eq. (|4.5| ) remains valid. In this way, even a CDM model can 
lead to an integrated Sachs Wolfe effect which then is known under the name 'Rees Sciama 
effect'. Furthermore, do to ultra violet radiation of the first objects formed by gravitational 
collapse, the universe might become reionized and electrons and radiation become coupled 
again. If this reionization happens early enough [z > 30) the subsequent collisions lead to 
additional damping of anisotropies on angular scales up to about 5°. However, present CMB 
anisotropy measurements do not support early reionization and the Rees Sciama effect is 
probably very small. Apart from these effects due to non-linearities in the matter distribu- 
tion, which depend on the details of the structure formation process, CMB anisotropies can 
be determined within linear perturbation theory. 

This is one of the main reason, why observations of CMB anisotropies may provide 
detailed information about the cosmological parameters (see contribution by A. Lasenby): 
The main physics is linear and well known and the anisotropies can thus be calculated within 
an accuracy of 1% or so. The detailed results do depend in several ways on the parameters 
of the cosmological model which can thus be determined by comparing calculations with 
observations. 

There is however one caveat: If the perturbations are induced by seeds (e.g. topologi- 
cal defects), the evolution of the seeds themselves is in general non-linear and complicated. 
Therefore, much less accurate predictions have been made so far for models where pertur- 



bations are induced by seeds (see, e.g., [[19], EO, |2Tj). In this case, the observation of CMB 



anisotropies might not help very much to constrain cosmological parameters, but it might 
contain very interesting information about the seeds, which according to present understand- 
ing originate from very high temperatures, T ~ 10 16 GeV. The CMB anisotropies might thus 
bury some 'fossils' of the very early universe, of the physics at an energy scale which we can 
never probe directly by accelerator experiments. 
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A An equation of motion for the Weyl tensor 



The Weyl tensor of a spacetime (M.,g) is defined by 

C^ p = i2% - 2«] + \ R 9 [ V ] P] , (Al) 

where [//... i/] denotes anti-symmetrization in the indices /i and v. The Weyl curvature has 
the same symmetries as the Riemann curvature and it is traceless. In addition the Weyl 
tensor is invariant under conformal transformations: 

C%( ff ) = C%,(a 2 <7) 

(Careful: This equation only holds for the given index position.) In four dimensional space- 
time, the Bianchi identities together with Einstein's equations yield equations of motion for 
the Weyl curvature. In four dimensions, the Bianchi identities, 

R/j,v[ap;\] = 

are equivalent to || 

C a ^ s ; s = W [a ^ ] - -g^ a R® . (A2) 
6 

This together with Einstein's equations yields 

C a ^ 5 ; s = 8vrG'(T 7[Q;/3] - -g^T^) , (A3) 

3 

where T pv is the energy momentum tensor, T = T x . 

Let us now choose some time-like unit vector field u, u 2 = —1. We then can decompose 
any tensor field into longitudinal and transverse components with respect to u. We define 

h\ = Q* v + u»u v , 

the projection onto the subspace of tangent space normal to u. The decomposition of the 
Weyl tensor yields its electric and magnetic contributions: 

= C pXua u x u a (A4) 
= \c, x ^u x rfiX ; (A5) 

where r] al3 ^ 5 denotes the totally antisymmetric 4 tensor with ??oi23 = \f~9- Due to symmetry 
properties and the tracelessness of the Weyl curvature, £ and B are symmetric and traceless, 
and they fully determine the Weyl curvature. One easily checks that £ pv and B pu are also 
conformally invariant. We now want to perform the corresponding decomposition for the 
energy momentum tensor of some arbitrary type of seed, T^ v . We define 



p s = T$vh? (A6) 

p s = l ~T$W» (A7) 

% = ~h;T^u a <H = -- a T m (A8) 

V = K h Mf-K uPs . (A9) 
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We then can write 

t Im) = Psu^Uy + psh^ + q^u u + ti M g„ + . (A10) 

This is the most general decomposition of a symmetric second rank tensor. It is usually 
interpreted as the energy momentum tensor of an imperfect fluid. In the frame of an observer 
moving with four velocity u, ps is the energy density, p$ is the isotropic pressure, q is the 
energy flux, u ■ q — 0, and r is the tensor of anisotropic stresses, T iW h ilu = r^u^ = 0. 

We now want to focus on a perturbed Friedmann universe. We therefore consider a four 
velocity field u which deviates only in first order from the Hubble flow: u = (l/a)<9o+ first 
order. Friedmann universes are conformally flat, and we require the seed to represent a 
small perturbation on a universe dominated by radiation and cold dark matter (CDM). The 
seed energy momentum tensor and the Weyl tensor are of thus of first order, and (up to 
first order) their decomposition does not depend on the choice of the first order contribution 
to u, they are gauge-invariant. But the decomposition of the dark matter depends on this 
choice. Cold dark matter is a pressure-less perfect fluid We can thus choose u to denote the 
energy flux of the dark matter, T^u u = —pcu^. Then the energy momentum tensor of the 
dark matter has the simple decomposition 

T^= Pc u,u v . (All) 

With this choice, the Einstein equations Eq. (|A3| ) linearized about an Q = 1 Friedmann 



background yield the following 'Maxwell equations' for E and B |22 

i) Constraint equations 



d'Bij = AirGrup^q^ (A12) 

-a p c D,j +-a p s ,j --a - - 
o 6 z a, 



d l £ t] = MG{\a 2 p c D, ] +\a 2 p S)3 -\dW ij -^q J ) . (A13) 



ii) Evolution equations 



a% + d%-a%% w ?/£ a ^ = -47TGa 2 h a{l r] jW ^T a ^ (A14) 
Sij + ^Sij + ah {i a rj j)MS u p B a r ' 5 = -^ixGiaqij - ^T {j + + apcUy), (A15) 



where denotes symmetrization in the indices % and j. The symmetric traceless tensor 

fields q^v and u^ u are defined by 



1 \ 

Q/J.V — ^nul -X 

1 x 

u fiu — u {n;v) T^h'fiuU -x • 



In Eqs. ( |A14| ) and ( A15 ) we have also used that for the dark matter perturbations only scalar 



perturbations are relevant, vector perturbations decay quickly. Therefore u is a gradient 
field, Ui = U-i for some suitably chosen function U. Hence the vorticity of the vector field u 
vanishes, = 0. With 

4 _ -2rpS A - -lrpS 

VOijk — O £-ijk , PS — a J-oo ana Qi — ~ a J-Qi ■, 
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we obtain from Eq. ( |A13| ) 



X 11" 



(A16) 



In Eq. ( |A16| ) and the following equations summation over double indices is understood, 
irrespective of their position. 

To obtain the equation of motion for the magnetic part of the Weyl curvature we take 
the time derivative of Eq. ( |A14|) , using u = (l/a)do + 1. order and r/oyfc = a 4 £ijk- This leads 
to 

a _ , „ r . a 



aft 



—a>{£lm(i[£j)l H £j)l)->m —4nG€i m (i[Tj)i, m H — Tj)i, m ]) 



(A17) 



a a 
where we have again used that u is a gradient field and thus terms like eijkUij,k vanish. We 
now insert Eq. (|A15 ) into the first square bracket above and replace product expressions of 
the form e^e^m and e^fcC/mn with double and triple Kronecker deltas. Finally we replace 
divergences of B with the help of Eq. (|A12 ). After some algebra, one obtains 



eim(i[£j)i + 



-V 2 B i:j - 4:7iGei m{i [2aqi 



a 
a 



Inserting this into Eq. (|A17 ) and using energy momentum conservation of the seed, we finally 
find the equation of motion for B: 



with 



q(B) r rpS j_ 

<~>ij — e lm(i[ I oiv)m' T : 



j)hm . 



(A18) 



(A19) 



Eq. ( |A18[ ) is the linearized wave equation for the magnetic part of the Weyl tensor in an 
expanding universe. A similar equation can also be derived for 8. 

Since dark matter just induces scalar perturbations and Bij is sourced by vector and 
tensor perturbations only, it is independent of the dark matter fluctuations. Equations 
Eqs. ( [AlCjj ) and ( |A18| ) connect the source terms in the Liouville equation of section 2, d' l £ i j 
and Bij to the perturbations of the energy momentum tensor. 
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